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Table 1. Parameters of the various collective modes 
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Table 11. Screening lengths A in directions perpendicular and parallel to the 
order parameter direction (z). 
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Table III. Coherence lengths characterizing different collective modes I 
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It was recently shown that conventional phonon-electron interactions may induce a triplet pairing 
state in time-reversal invariant 3D Dirac semi - metals. Starting from the microscopic model of 
the isotropic Dirac semi-metal, the Ginzburg-Landau equations for the vector order parameter is 
derived using the Gor’kov technique. The collective modes including gapless Goldstone modes, and 
gapped Higgs modes of various polarizations are identihed. They are somewhat analogous to the 
modes in the B phase of He^, although in the present case quantitatively there is a pronouneced 
difference between longitudinal and transverse components. The difference is caused by the vector 
nature of the order parameter leading to two different coherence lengths or penetration depths. 
The system is predicted to be highly dissipative due to the Goldstone modes. The time dependent 
Ginzburg - Landau model in the presence of external fields is used to investigate some optical 
and magnetic properties of such superconductors. The AG conductivity of a clean sample depends 
on the orientation of the order parameter. It is demonstrated that the difference between the 
penetration depths results in rotation of the polarization vector of microwave passing a slab made 
of this material. The upper critical magnetic field 11^2 was found. It turns out that at fields close 
to Hc 2 the order parameter orients itself perpendicular to the field direction. In certain range of 
parameters the triplet superconducting phase persists at arbitrarily high magnetic field like in some 
p wave superconductors. 


I. INTRODUCTION 

Recently 3 D Dirac semi - metals (DSM) like Na^Bi and Cd^As2 with electronic states described by Bloch wave 
functions, obe;^ng the ” pseudo-relativistic” Dirac equation (with the Fermi velocity vp replacing the velocity of light) 
were observed^ and attracted widespread attention. The discovery of the 3 D Dirac materials makes it possible to 
study their physics including remarkable electronic properties. This is rich in new phenomena like giant diamagnetism 
that diverges logarithmically when the chemical potential approaches the 3 D Dirac point, a linear-in-frequency AC 
conductivity that has an imaginary partQ, quantum magnetoresistance showing linear field dependence in the bulkQ. 
Most of the properties of these new materials were measured at relatively high temperatures. However recent exper¬ 
iments at low temperature on topological insulators and suspected 3 D Dirac semi-metals exhibit superconductivity. 
Early attempts to either induce or discover superconductivity in Dirac materials were promising. The well known 
topological insulator 312863 doped with Cu, becomes superconducting at Tc = 3 . 8 iT[ 5 . At present its pairing sym¬ 
metry is unknown. Some experimental evidence^ point to a conventional phononic pairing mechanism. The spin 
independent part of the effective electron - electron interaction due to phonons was studied theoreticallyFor a 
conventional parabolic dispersion relation, typically independent of spin, the phonon mechanism leads to the s-wave 
superconductivity. The layered, non-centrosymmetric heavy element PbTaS62 was found to be superconducting [^. 
Its electronic properties like specific heat, electrical resistivity, and magnetic-susceptibility indicate that PbTaS62 is 
a moderately coupled, type-II BCS superconductor with large electron-phonon coupling constant of A = 0 . 74 . It was 
shown theoretically to possess a very asymmetric 3 D Dirac point created by strong spin-orbit coupling. If the 3 D is 
confirmed, it might indicate that the superconductivity is conventional phonon mediated. 

More recently when the Cu doped Bi2Se3 was subjected to pressureRl, Tc increased to 7 K at SOGPa. Quasilinear 
temperature dependence of the upper critical field Hc2, exceeding the orbital and Pauli limits for the singlet pairing, 
points to the triplet superconductivity. The band structure of the superconducting compounds is apparently not very 
different from its parent compound 312863, so that one can keep the two band k • p description (S'e Pz orbitals on the 
top and bottom layer of the unit cell mixed with its neighboring Bi pz orbital). Electronic-structure calculations and 
experiments on the compounds under pressure^ reveal a single bulk three-dimensional Dirac cone like in Bi with 
large spin-orbit coupling. Moreover very recently some pnictides were identified as exhibiting Dirac spectrum. This 
effort recently culminated in discovery of superconductivity in C(i3As2 0 . It is claimed that the superconductivity is 
p-wave at least on the surface. 

The case of the Dirac semi-metals is very special due to the strong spin dependence of the itinerant electrons’ effective 
Hamiltonian. It was pointed out@ that in this case the triplet possibility can arise although the triplet gap is smaller 
than that of the singlet, the difference sometimes is not large for spin independent electron - electron interactions. 
Very recently the spin dependent part of the phonon induced electron - electron interaction was considered (l^ and it 
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was shown that the singlet is still favored over the triplet pairing. Another essential spin dependent effective electron- 
electron interaction is the Stoner exchange among itinerant electrons leading to ferromagnetism in transition metals. 
While in the best 3D Weyl semi-metal candidates it is too small to form a ferromagnetic state, it might be important 
to determine the nature of the superconducting condensate. It turns out that it favors the triplet pairing[ll|. Also a 
modest concentration of magnetic impurities makes the triplet ground state stable. 

In a multicomponent superconductor collective modes including gapless Goldstone modes (sound), and gapped 
Higgs modes of various polarizations play an important role in determining thermal and optical properties of the 
materialfl^. In addition, as mentioned above, generally the applied magnetic field is an ultimate technique to probe 
the superconducting state. In a growing number of experiments, in addition to magnetotransport, magnetization 
curves, the magnetic penetration depth and upper critical ma gnet ic field were measured[l3. It is therefore of im¬ 
portance to construct a Ginzburg - Landau (GL) description [ij of these novel materials. This allows to study 
inhomogeneous order parameter configurations (junctions, boundaries, etc.), the collective modes (somewhat analo¬ 
gous to the modes in the B phase of He^ and magnetic and optical response that typically involve inhomogeneous 
configurations (like vortices) not amenable to a microscopic description. 

In the present paper we derive such a GL type theory for triplet superconductor from the microscopic isotropic 
DSM model with attractive local interaction. The order parameter in this case is a vector field and the GL theory of 
vector field already considered in literaturein connection with putative p - wave superconductors have several 
extraordinary features, both quantitative and qualitative. 

The paper is organized as follows. The model of the (phonon mediated or unconventional) local interactions of 3D 
Dirac fermion is presented and the method of its solution (in the Gorkov equations form) including the symmetry 
analysis of possible pairing channels and the vectorial nature of the triplet order parameter is given in Section 11. 
In Section III the Gorkov formalism, sufficiently general to derive the GL equations, is briefly presented. The most 
general form of the GL energy of the triplet superconductor in magnetic field consistent with the symmetries is given 
in IV. The coefficient of the relevant terms are calculated from the microscopic DSM model in section V. Section 
VI is devoted to applications of the GL model. The ground state degeneracy, the character of its excitations and 
basic magnetic properties are discussed. The vector order parameter is akin to optical phonons with sharp distinction 
between transverse and longitudinal modes. Transverse and longitudinal coherence lengths and penetration depths 
are calculated and the upper critical magnetic field is discussed. Section VI includes generalizations to include Pauli 
paramagnetism, discussion of an experimental possibility of observation of the excitation and conclusion. 


II. THE LOCAL PAIRING MODEL IN THE DIRAC SEMI-METAL. 


A. Pairing Hamiltonian in the Dirac semi-metal. 


Electrons in the 3D Dirac semimetal are described by field operators ipfs (r), where f = L,R are the valley index 
(pseudospin) for the left/right chirality bands with spin projections taking the values s =t 7 i with respect to, for 
example, the 2 axis. To use the Dirac (’’pseudo-relativistic”) notations, these are combined into a four component bi¬ 
spinor creation operator, '0^ = ^ whose index 7 = {/, s} takes four values. The non-interacting 

massless Hamiltonian with Fermi velocity vp and chemical potential fi reads [l^ 


K 


(r) K-ip (r); 


-ihvFV''al^g - fi6^s, 


where the three 4x4 matrices, i = x,y,z, 


ot = 




( 1 ) 

( 2 ) 

( 3 ) 


are presented in the block form via Pauli matrices cr. They are related to the Dirac 7 matrices (in the chiral 
representation, sometimes termed ’’spinor”) by a = P'y with 


/3 = 


Here 1 is the 2x2 identity matrix. 


0 1 
1 0 


( 4 ) 
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We consider a special case of 3D rotational symmetry that in particular has an isotropic Fermi velocity. Moreover 
we assume time reversal, Qtp (r) = iaytp* (r), and inversion symmetries although the pseudo- Lorentz symmetry will 
be explicitly broken by interactions. The spectrum of single particle excitations is linear. The chemical potential fj, is 
counted from the Dirac point. 

As usual in certain cases the actual interaction can be approximated by a model local one: 


^eff = - 


9 


(r) (r) (r) (r). 


Unlike the free Hamiltonian K, Eq. o, this interaction Hamiltonian does not mix different spin components. 
Spin density in Dirac semi-metal has the form 

S (r) = (r) (r), 


where the matrices 


s= -OL 75 = 

75 = 


CT 0 
0 (T 


-1 0 

0 1 


are also the rotation generators. 


( 5 ) 


( 6 ) 

( 7 ) 


B. The symmetry classification of possible pairing channels. 


Since we consider the local interactions as dominant, the superconducting condensate (the off-diagonal order pa¬ 
rameter) will be local 

O = ^ V'i (r) (r), (8) 

where the constant matrix M should be a 4 x 4 antisymmetric matrix. Due to the rotation symmetry they transform 
covariantly under infinitesimal rotations generated by the spin S'® operator, Eq.ljH): 


f 4 ’t (r) (r), (r) {r') 

J v.r' ^ 

- / V-+ (r) (r). 


(9) 


Here and in what follows denotes the transpose matrix. The representations of the rotation group therefore 
characterize various possible superconducting phases. 

Out of 16 matrices of the four dimensional Clifford algebra six are antisymmetric and one finds one vector and 
three scalar multiplets of the rotation group. The multiplets contain: 

(i) a triplet of order parameters: 


The algebra is 


= {-P^z, -iplz, /3ax} = {Tx, Ty, T^} 
^iTj + TjTj\ = 2i£ijkTk. 


Note that the three matrices Ti are Hermitian. 
(ii) three singlets 


( 10 ) 


( 11 ) 


M^ = iay] M^ = iEy; =-iPttyj^. ( 12 ) 

Which one of the condensates is realized at zero temperature is determined by the parameters of the Hamiltonian and 
is addressed next within the Gaussian approximation. As was shown in our previous work [Til . El, either exchange 
interactions or magnetic impurities make the triplet state a leading superconducting channel in these materials. 
Therefore we will consider in the next section only the vector channel. 
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III. GORKOV EQUATIONS AND THE TRIPLET PAIRING 
A. Gorkov equations for Green’s functions in matrix form 

Using the standard BCS formalism, the Matsubara Green’s functions (r is the Matsubara time) 

Ga/3 (r,T;r',r') = - (r, r) V’js (r', t')) ; 


(13) 


obey the Gor’kov equations [^: 


^b /3 (r,T;r',T') = ^ {r, t) ^^ {r', t')) ; 


dG^K (r,r;r',T') 
dr 




L ^ 

Kjp 


(r",r;r',r') 
(r — r') 6 {t — t' 


(14) 


-5^/37 (J*, r- r, r) F+^ (r, r, r', r') = (r - r') d (r - r') ; 
9F+^ (r,T;r',T') 


(9r 

(r, r, t) G/3k (r, t, r' , t') = 0. 

These equations are conveniently presented in matrix form (superscript t denotes transposed and / - the identity 
matrix): 


D-i {X,X”)G{X”,X')- 
-A{X)F+ {X,X') 


(15) 


= I5{X-X')- 


[ (A, A") F+ (A", A') + A** (A) G (A, A') = 0. 

Jx" 


Here A = (r,r), (A) = gFp^ {X,X) and 


d 


D-;{X,X') = -5^p—5{X-X')- 


(16) 


-(5(r-r')(r 


Fa.p 


In the homogeneous case the Gor’kov equations for Fourier components of the Green’s functions simplify consider¬ 
ably: 


D ^ {u},p)G{uj,p)-AF+{u},p) = /; (17) 

(a;,p)F+(a;,p)-k A**G(a;,p) = 0. 

The matrix gap function can be chosen as 

Xp'j = gF^p (0) = A^M^p, (18) 

with real constant A^. Here D~^ {uj, p) = -I- /i — a • p, is the noninteracting inverse Dirac Green’s function for the 

Hamiltonian Eq. © and D ^ (w, p) = iu) — g, — a*-p, where = ttT (2n + 1) is the fermionic Matsubara frequency. 
Solving these equations one obtains (in matrix form) 

G-i = D-^+ADA**; (19) 

F+ = -5a**g, 

with the gap function to be found from the consistency condition 

A‘* = -g^bX*G. 

ujp 

Now we find solutions of this equation for each of the possible superconducting phases. 


( 20 ) 
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B. Homogeneous triplet solution of the gap equation. 

In this phase rotational symmetry is spontaneously broken simultaneously with the electric charge U (1) (global 
gauge invariance) symmetry. Assuming 2 direction of the p - wave condensate the order parameter matrix takes a 
form: 


A = 


( 21 ) 


where A^ is a constant. The energy scale will be set by the Debye cutoff Tp of the electron - phonon interactions, 
see below. 

The spectrum of elementary excitations at zero temperature was discussed in ref. [lll |. There is a saddle point 
with energy gap 2A2 on the circle + Py = p'^/vp, Pz = 0. The gap A^ as a function of the dimensionless 

phonon-electron coupling A = gN, where N being the density of states (all spins and valleys included), increases upon 
reduction in p. At large p » Tp, as in BCS, the gap becomes independent of p and one has the relation 


1 

g 


-sinh-i 

12 



2 p^ 

Tr’^VphP ’ 


( 22 ) 


leading to an exponential gap dependence on A when it is small: Aj, = Tp/sinh (12/A) ~ 2Tp,e 

The critical temperature is obtained from Eg . (1201) with disctret uj by substituting A^ = 0. To utilize the orthonor¬ 
mality of Ti, Ti[TiT*] =4%, one multiplies the gap equation by the matrix T^jg and takes the trace: 


The bubble integral is 


^Tr (T,r/) = - = TJ5z. 
9 9 




^ Tr (t.DT/d) = 4%Te x 


E 

np 


(pi - pI) + + 


+ 2^2 (vj,p2 + ^ 2 ) 


(23) 


(24) 


Performing first the sum over Matsubara frequencies and then integrate over q one obtains, similarly to the singlet 
BCS, (see Appendix A for details): 


r,= 




(25) 


where logy^; = 0.577 is the Euler constant. 


IV. A GENERAL GL DESCRIPTION OF A TRIPLET SUPERCONDUCTOR IN A MAGNETIC FIELD. 

In this section the effective description of the superconducting condensate in terms of the varying (on the mesoscopic 
scale) order complex parameter vector field A^ (r) is presented. 


A. The GL description for a vector order parameter 

The static phenomenological description is determined by the GL free energy functional F [A (r), A (r)] expanded 
to second order in gradients and fourth order in A. In a magnetic field B, as usual, space derivatives of the 
microscopic Hamiltonian become covariant derivatives V — e* = 2e due to gauge invariance under 
Ai —>■ Ai, Ai —>■ Ai — IfVy. Naively the only modification of the GL energy is in the gradient term, Ea. d^Sl) : 

the most general gradient term consistent with rotation symmetry and the U (1) gauge symmetry is 


grad 



Ut 


nVyA,nv,A,)-\ ] 
\-(AA,)*(D,A,) / \ 
+UL (AAj)* {V,A,) J 


(26) 
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The factor N, the density of states, is customarily introduced into energy [HI- It was noted in [l^ . that, unlike in 
the usual scalar order parameter case, the longitudinal and transverse coefficients are in general different, leading to 
two distinct coherence lengths, see Section V. Possibilities for the local terms are[l5l| 

Floe = |a (T - To) A*A, + ^ (A*A,)' + ^ |A,A,|"|. (27) 

The magnetic part, F^ag = completes the free energy. 

B. The set of time independent GL eqnations for triplet order parameter 

The set of the GL equations corresponding to this energy are obtained by variation with respect to 'k* and At. The 
first is: 


- |wT i + \nL A,+ (28) 

+a (T - To) A, + /3iA,A*Aj + = 0. 

The anticommutator appears due complex conjugate terms in Eq. (iMiidll- The Maxwell equation for the supercurrent 
density is: 


ie* 

J^ = —N {utA*P,Aj + uA*VjA,) + cc, (29) 

where u = ul — ut- 

Having defined coefficients ut,l,I3i 2 a, our aim in the next Section is to deduce them from the microscopic 
Dirac semi-metal model. 


V. GL COEFFICIENTS FROM THE GOR’KOV EQUATIONS 

For the calculation of coefficients of the local part, the homogeneous Gor’kov equation, Ea. (l20l) suffices, while for 
calculation of the gradient terms a general linearized equation, Eq. (USD is necessary. 


A. Local (potential) terms in Gor’kov 

Iterating once the equation Eq.([20D with help of Ea. dT^ one obtains the local terms to third order in the gap 
function: 


1 A** 
9 

+E 

Ljp 


(30) 


D{cu,p)A**D^{io,p)- 
—D {uj, p) A**D {uj, p) AD {oj, p) A**D (w, p) 


Using A** = A*Ti, multiplying by T/ and taking the trace, one gets the linear local terms 


Na {T-To)A* = -A*, 


(31) 


where the bubble integral was given in Ea. (l24D . Expressing g via To, see Ea. (l25ll allows to write the coefficient a of 
A* in the Gorkov equation Ea. d^ is 


irr. 8 m" , T AT-To 

a{T- To) = log ; 


^To i To 


(32) 
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The cubic terms in Ea. (l30)) . multiplied again by T/ and ’’traced” take the form 

N{f3,A*A,A*+f],A*A*A,) 


= -A*AfcAr^Tr{r/ 

cjp 

The calculation is given in Appendix A and results in: 

7C(3) 1 . 


'IDT" DTkDTl* D 


}■ 


Pi = 


207r2 T2’ 


(33) 


(34) 


The Riemann zeta function is C (3) = 1.2. 


B. Linear gradient terms 

To calculate the gradient terms, one first linearizes the Gor’kov equations, Ea. dTSl) 


f D-i(A,A")G(A",A') 

Jx" 

= IS{X - X') ^ G = D-^; 

F+ (A, A') 

= - [ £)* (A - A") A*‘(A") (A" - A'). 

Jx" 


(35) 


In particular, 


^A** (A) 
9 


F+ (A, A) 



A') A** (A') D (A' - A). 


(36) 


The anomalous Green’s functions are no longer space translation invariant, so that the following Fourier trans¬ 
form is required: The (time independent) order parameter is also represented via Fourier components A* (A) = 
e“®^ '’A* (P). The linear part Gor’kov equation (this time including nonlocal parts) depending on the ’’external” 
momentum P reads: 


iA*‘ {V) + Y,D (u;,p) A‘* (P) G (a;,p - P). (37) 

9 

To find the coefficients of the gradient terms, one should consider contributions quadratic in P from the expansion 
of both A** (P) and H (w,p — P). In view of the gap equation Ea. (l23l24l) . the expansions of A‘* (P) cancel each 
other up to small corrections of order T — T^. So that multiplying by Tf and taking the trace 

\PuPi Y. Tr {t/5 { uj , p) Tf 74" (w, p)} A*, (38) 


where D'h (w,p) = Comparing this with the gradient terms in the GL equation, Ea. (l28l) . see Appendix B 

for details, one deduces 


Ut 


28C{i)vl,h^^ _ 1 

157r2 Tf ’ X 


(39) 


Note the very small longitudinal coefficient, ul « ut- As we shall see in the following section it has profound 
phenomenological consequences. 







VI. BASIC PROPERTIES OF THE TRIPLET SUPERCONDUCTOR 


A. Ground state structure and degeneracy 

A ground state is characterized by three independent parameters corresponding to three Goldstone bosons. The GL 
energy is invariant under both the vector 0(3) space rotations, —>■ RijAj, and the superconducting phase U{1), 
Ai —In the superconducting state characterized by the vector order parameter A (|A| = A, energy gap) 
the U (1) is broken: U (1) —>■ 1, while the 0(3) is only partially broken down to its 0(2). There are therefore three 
Goldstone modes. Here we explicitly parametrize these degrees of freedom by phases following ref.jT^. Generally a 
complex vector field can be written as 


A = A (n cos x+im sin y), (40) 

where n and m are arbitrary unit vectors and 0 < x < 7r/2, see Fig. 1. 

Using this parametrization the homogeneous part of the free-energy density, Eg. (ITT)) , takes the form 

^ f a(T-r,)A“ + iftA-+ 

JV I (cos^ (2x) + (n • m)^ sin^ (2x) j 

This form allows us to make several interesting observations. The crucial sign is that of /32- In previous studies [Tsll^ 
only /?2 > 0 (so called phase A) was considered. In our case however /32 < 0 and different ground state configurations 
should be considered. In phase B the minimization gives, n = ±m. Note two different solutions. So that the ’’vacuum 
manifold” is 



A = Aone*x. (42) 

Here the range of x was enlarged, —7r/2 < x < to incorporate n = ±m. The ground state energy density 
therefore is achieved at 


2 _ Q (Te - T) _ a (Tg - T) 

° /3i+/32 /3 ■ ^ ^ 

Mathematically the vacuum manifold in phase B is isomorphic to S 2 ® S\lZ 2 .. This determines the thermodynamics 
of the superconductor very much in analogy with the scalar superconductor with B = Pi+ j32- However the collective 
modes, the AC conductivity and the magnetic properties are markedly different [23|. 


B. Collective excitation modes. 

Here the response of the superconductor in phase B to an external perturbation, like boundary or magnetic field, 
is considered. The basic excitation modes are uncovered by the linear stability analysis very similar to the so-called 
Anderson - Higgs mechanism in field theory applied to (scalar order parameter) superconductivity a long time ago 
[ 2 ^. Two basic scales, the coherence length (scale of variations of the order parameter) and magnetic penetration 
depth (scale of variations of the magnetic field), are obtained from the expansion of the GL energy to second order 
in fluctuations around superconducting ground state at zero field. In the superconducting state the most convenient 
gauge is the ’’unitary” gauge in which the U (1) phase of the order parameter is set to zero, so that we are left with 
three massive vector potential fields Ai,A 2 ^A^. Ghosing the grond state as A = Aq (0,0,1), see Fig.2 the order 
parameter in the unitary gauge can be parametrized by five real fields 


A = Ao(I+e)(i?i+iJi,i?2+i/2,l). (44) 

The sponteneous breaking of the space rotation symmetry O (3) into its O (2) subgroup of rotations within the x — y 
plane according to the Goldstone theorem leads to two gapless modes JTq,, a = 1,2. Due to the residual symmetry the 
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Fourier components of the fluctuation fields can be generally written as combination of the radial and the tangential 
components: 


Ra — (^Rr^a 

Ice — {Irka “f It^ajskjS^ /^_L 


e u, 


1/2 


he 


{A^kcK ^t^ajsk^^ (k_\_j 


as = 


e u 


he 


1/2 


(45) 

(46) 


(47) 


1 /2 

with k± = (fc^ + ky) . The GL energy to quadratic order in eight fluctuation fields rj = {e, R, It, Rr, Rt, ar, at, 03 }, 




(48) 


decomposes into three independent sectors, where M is dimensionless fluctuation matrix. 

1. Three massive fields R, ara^ can mix due to kinetic terms: 


Ml 


/ 

V 


2^0 1 /^ 2 ! + 

+UTk'^ + ukj_ 

i^ks 

Urji 

0 


-i^ks 0 \ 

Urji 

1 “h X'j'k‘^ — X'j'k^k±^ 

— X'j'k^k±^ URj X'Y'k‘^ / 


where /S7ruTe'^‘^NA q. 

To order the eigenvalues are, 


E {k^) jv = (n^ + dkl + . 


(49) 


(50) 


The values of the gaps in the spectrum 17^ and and corresponding velocities in directions perpendicular and parallel 
to the vector order parameter (taken to be z), C±, Cy are given in Table I 


2. The tangential massive, It, at, sector 


/ 2 I/32I 2 ^ 0 i-^ks \ 

3 . The Goldstone transverse mode Rt does not mix, M 3 = urk"^, so that Eq/V = AAqWtA:^. 

4. Higgs and Goldstone radial collective modes £, Rr form a 2 x 2 matrix: 

,, _ / urk"^ + uk\ uk±k3 \ 

^ \ uk±k3 urk'^+uk^ + Ad ) 


(51) 


(52) 


C. Coherence length and penetration depth for massive collective modes 

The six ’’massive” fields, e,Ia and A with different longitudinal and transversal characteristic lengths: l± = 
l|| = C\\/Il ,and the same for ^ and A, see Table 1. 

This is different compared with the one component (singlet) superconductor in two respects. First the number of 
Higgs modes is larger since in addition to the superfluid density determined by s, there are two additional components 
R and R. Second, as mentioned above, since the superconducting condensate is oriented, one has two different 
velocities. As far as (massive) photon modes are concerned, the number of modes remains the same but the anisotropy 
persists. Let us start with m 
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Here definitions of the coherence lengths (of fluctuation of the superfluid density e perpendicular and parallel to 
the direction of the order parameter n), correlation lengths of the relative weights between different components of 
the vector order parameter {Ir,It) and screening lengths (vector potential A fluctuations): 


\l 


‘2ul f.2 _ Ut ^2 
UT ,2 _'^L, 

til - -t 


2 1/321 Ag’ 


STruTe*^ AAg 


UL ,2 . 
Ut 


,2 _ Ut ,2 

Al — - Arp. 

UL 


(53) 


In the isotropic superconductor one recovers the standard formulas since ut = ul 

Our calculation in the previous Section for the Dirac semi-metal, see Eg. (1391) . demonstrate that both are quite 
different since ul/ut = 1/32 << 1. This is obviously of great importance for large magnetic field properties of such 
superconductors and will be discussed below. 


VII. MAGNETIC AND OPTICAL PROPERTIES 


A. Strong magnetic fields: is there an upper critical field Rc2? 


In strong homogeneous magnetic field H (assumed to be directed along the z axis) superconductivity typically (but 
not always, see an example of the p-wave superconductor that develops flux phases (23|) disappears at certain critical 
value Hc 2 - This bifurcation point is determined within the GL framework by the lowest eigenvalue of the linearized 
GL equations. This is an exact requirement of stability of the normal phasefl^, The linearized GL equation 
Eg. d^ reads: 


(a - utV^) Sij 


■ {Vi,Vj} 


A, = 0, 


(54) 


where coefficients are in Eo. (l39l) . and u = ul — ut- We use the Landau gauge, = Hc 2 y', Ay = Az = 0. Assuming 
translation symmetry along the field direction, dzAi = 0, the operators of the eigenvalue problem depend on x and y 
only. 

Since we have three components of the order parameter, there are three eigenvalues. It is easily seen from Eg. (1541) 
that the z- component of the order parameter parallel to the external field direction is independent of the other 
two, Aj,, Ay, leading to the ordinary Abrikosov value: 


- utV'^Az = -aAz -t = 


^0 

27r^y ’ 


(55) 


where = ut (see Eg. (1391) . To avoid confusion with customary notations for layered materials (like high Tc 
cuprates), the material that is modelled here is isotropic and ’’parallel”, ’’perpendicular” and refer to the relative 
orientation of the magnetic field to the vector order parameter rather than to a layer. The orientation of the order 
parameter in isotropic material considered here, due to degeneracy of the ground state, is determined by the external 
magnetic field as we exemplify next. 

The two remaining eigenvalues involving only the order parameter components A^, and A^ perpendicular to the 
field (see Fig.3) are obtained from diagonalizing the ’’Hamiltonian”: 


n 


A, 

A. 


A,, 

A„ 


_ _ / utVI + ulVI 


(56) 


utVI -I- ulVI 


This nontrivial eigenvalue problem fortunately can be solved exactly, see Appendix C. The lowest eigenstate being 
a superposition of just two lowest even Landau levels, |0) and |2) are given. The lowest of these eigenvalues is 




I {ut + Ul) _ 

-•^3 (u|, + ul) - 2 utUl 


= a{Tc-T). 


he 


(57) 
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The corresponding critical field (’’perpendicular” refers to the order parameter direction (Fig.3) ) that can be 
expressed via an effective ” perpendicular” coherence length, 


H^2 = 


d>n 




(58) 


It is always larger than therefore is physically realized. The upper field becomes infinite at Tc = 

ul/ut = ( l 3 — 4 -\/ T 0 ) /3 ~ 0 . 117 . This means that in such material superconductivity persists at any magnetic field 
like in some p- wave superconductors. It was found in Section IV that for the simplest Dirac semi-metal, r = 1/32 < Tc 
see Eg. (15^ . Thus there is no upper critical field in this case. Of course, different microscopic models that belong 
to the same universality class, might have higher r. In any case the Abrikosov lattice is expected to be markedly 
different from the conventional one and even from the vector order parameter model studied in [T5j |. 


B. Dissipative dynamics 

The set of the GL equations corresponding to this energy are obtained by variation with respect to A* and Ai. The 
first is the time dependent GL equation (the covariant derivative is replaced by partial since in a superconductor the 
scalar potential can be taken to be zero on the mesoscopic scale): 


- ra* A, 



(59) 


This should be supplemented by the Maxwell equation including the normal metal contribution the the current 
J" = J — J'*, J® = (uTA*'DiAj + uA*'DjAi) + cc. This determines the dynamics of the vector potential: 


-^dtA, = (60) 

In the small fluctuations approximation the dominant role is played by the two Goldstone modes, Rt,Rr, due to 
spontaneous breaking of the 3D rotation symmetry. The Rt mode is still isotropic, while Rr is not Neglecting the 
massive excitations the dissipative dynamics of the diffusion type is governed by 

dtRt = Drk^Rt; (61) 

dtRr = {DTkj_ + DLkl) Rr, 


where Dt,l = The diffusion coefficient of this equation is anisotropic and is discussed in Section VIII. This 

would lead to increase in thermal conductivity inside the superconducting state even at low temperature. 


C. The AC conductivity 

In external AG field represented by (no spatial dispersion), A = (w) one obtains in linear response 


J^ = 


2ie*^NAl 

ojTP- 


{urEi + UE3631) + anEi. 


Therefore the conductivity tensor reads: 




CTn — Cr^ 0 0 

0 an — a?p 0 

0 0 an - a^ 


(62) 


(63) 


where a?p (w) = The AG conductivity this is different for the order parameter and the perpendicular 

directions. One of the interesting consequences of this phenomenon is rotation of the polarization of microwave that 
passes the DSM film. 
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D. Rotation of the polarization of the microwave 


The material becomes optically active, i,e. the polarization of the electromagnetic wave rotates. The dispersion 
relation is: 


Attlo 


(a;)B = B. 


(64) 


■ Z iZ 

The perpendicular to the order parameters are eigenvectors with eigenvalue -- ——j-j, while the third eigenvector 

• 2 1 2 

(0, 0,1) has the eigenvalue is ——j-j. Assume that the incident electromagnetic wave described on the surface 
of DSM by the magnetic field Bq is perpendicular to the order parameter direction taken as z (see Fig.4). 


Without loss of generality it can be taken as x (due to the residual 0(2) rotation symmetry). The Fourier 
component of the magnetic field in the y — z plane are 

By = Bo exp (65) 

Bz = Bo exp . 

The complex wave vectors are 







1 + 




~~2 


+ 1 ( \ \l\ 

~~2 


( 66 ) 


1/2 

Here the screening length is defined as usual J = (c/47ra;(T„) ' . 

The penetration of the microwave radiation into the sample surface demonstrates differences for AC component 
parallel and perpendicular to the direction of the order parameter resulting in effective rotation of the incident wave 
vector like at Faraday effect. 


VIII. DISCUSSION AND CONCLUSIONS 
A. The vector nature of the order parameter 

The physical properties of the triplet superconductivity appearing in 3D Dirac semi-metals were considered. Start¬ 
ing from the microscopic model of the isotropic Dirac semi-metal, the Ginzburg-Landau energy for this field is derived 
using the Gor’kov technique. The properties of the triplet superconductor phase of the Dirac semi-metal has ex¬ 
tremely unusual features that we would like to associate qualitatively with the characteristics of the Gooper pair. 
The superconducting state generally is a Bose - Einstein condensate of composite bosons - Gooper pairs, classically 
described by the Ginzburg - Landau energy as a functional of the order parameter. In the present case the Gooper 
boson is described by a vector field (r). In this respect it is reminiscent to phonon and vector mesons in particle 
physics [i^. 

Vector fields generally have both the orbital and internal degrees of freedom often called polarization. The internal 
degree of freedom might be connected to the ’’valley” degree of freedom of constituents of the composite boson. We 
have provided evidence that the Gooper pair in DSM has finite orbital momentum, albeit, as will be shown shortly, the 
spin magnetic moment is zero. Microscopically the unusual nature is related to the presence of the valley degeneracy 
in Dirac semi-metal. While in a single band superconductor the Pauli principle requires a triplet Gooper pair to have 
both odd angular momentum and spin, it is no longer the case in the Dirac semi-metal. 

A massive bosonic vector field in isotropic situation (the case considered here) generally have distinct transversal 
and longitudinal polarizations (massless fields like photons in dielectric do not possess the longitudinal degree of 
freedom). The results for collective modes in triplet superconductor in DSM demonstrate pronounced disparity 
between dispersion of various polarizations, see Table I. In particular we have found sound velocity of the two Goldstone 
modes and, screening lengths of three gapped photon modes and three coherence lengths of the other (Higgs) gapped 
modes. These all have an impact on transport, optical and magnetic properties of these superconductors. 
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B. Estimates of the characteristics of the collective modes and Faraday effect in a typical Dirac semi-metal 


Substituting the values of parameter of the vector GL equation found in Section IV into formulas for various 
coherence lengths described in Section VI, one obtains, 




7( (3) ^2 _ ooc2 

807r2 T2 ’ 


32^1; 




2 _ 
02 — 




37r 

e*2/i2 ’ 


^OL 


— 32Agp, 


(67) 


where (T) =Co_L/(l-r/rc) etc. 

For a typical DSM one estimates the Fermi velocity and chemical potential [i| = c/200, fj, = 0.2eV, and with 

the expected critical temperature Q, [I3| Tc = 5K, one obtains for = 230nm and Apr = 220nm. This has an 
impact on the magnetic flux penetration into this kind of superconductors. The value of Abrikosov parameter k = X/^ 
depends on mutual direction of the DC magnetic field. While kt = Xt/^± ~ 1, = Xp/^p ~ 30. 

Let us estimate the characteristics of the two Goldstone modes, arising in the triplet superconducting state due to 
spontaneous breaking of the rotational O (3) symmetry. The isotropic transverse mode Rt defined in Section III is 
isotropic and its dynamics is described by Ea. (l6ip with diffusion constant Dp = Nup/T, while the anisotropic mode 
Rr involves both the transverse and the longitudinal constant that is different. Dp = Nup/T. To estimate these, let 
us exploit the relation [l^ the time constant as F = Thus 


Dt 


Dl 


25 • 7C (3) vlh ^ 2 ^ 

IStT^ Tc s 

7((3)vpti 
IStT^ Tc s 


( 68 ) 

(69) 


These considerations were made for superconductor without significant pinning - disorder on the mesoscopic scale. 
Recently the AC response of the disordered superconductor was utilized to probe Goldstone modes [2^. We have 
demonstrated that they are abundant in the triplet DSM superconductor. Therefore it would be interesting to look 
for effects of the Goldstone modes including damping resulting in strong sound absorption in these systems. 

We have calculated the AC conductivity of DSM, see Ea. dMl) . and applied it to describe an intriguing effect of 
optical activity of DSM. According to Ea. (|66|) . the polarization vector of the incident beam rotates while passing a 
film of thickness d by (see Fig.4) 

tan^ = \Bz (d) /By (d)] = exp [(A^^ — A/^) d] , (70) 


1 /2 

under the assumption the skin depth 6 = (cjTKijjGr/) ' is much larger than both At and Xp. 
Xp = 32At, At = 220nTO,and d = 1mm, </>—>■ tt/ 2 significantly deviation from the initial 7r/4 value, 
the effect is due to the difference between the two penetration depths. 


In particular for 
Let us stress that 


C. DSM superconductor under constant magnetic field 

Several new features appear when an external field is applied. The Ginzburg - Landau model was used to determine 
upper magnetic field Hc 2 - It turns out that the lowest energy solution is when the order parameter of the texture 
orients itself perpendicular to the field direction. We have shown that the upper field becomes infinite for up/up < 
(13 — 4-\/l()) /3 ~ 0.117. In particular it is obeyed within our microscopic model. This means that in such material 
superconductivity persists at any magnetic field like in some p- wave superconductors. The expression for the Hc 2 for 
up/uT > 0.117 is given by Ea. (|^ . 

The vortex physics of strongly type II triplet superconductors of this type is very rich and some of it has already 
been investigated in connection with heavy fermion and other superconductors suspected to possess p-wave pairing. 
In particular, their magnetic vortices appear as either vector vortices or so-called skvrmions[l^ - coreless topologically 
nontrivial textures. The magnetic properties like the magnetization are very peculiar and even without a magnetic field 
the system forms a ’’spontaneous flux state”. The material therefore can be called a ’’ferromagnetic superconductor”. 
The superconducting state develops weak ferromagnetism and a system of alternating magnetic domains It was 
noted [l^ that the phase is reminiscent to the phase B of superfluid He^. [3 (with an obvious distinction that the 
order parameter in the later case is neutral rather than charged and tensorial rather than vectorial). 
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Since the prediction of the FFLO effect [T^ in low Tc superconductors it is well known that at very high magnetic 
fields the direct spin - magnetic field coupling on the microscopic level might not be negligible. The singlet channel 
Cooper pair is effectively ’’broken” by the splitting since the spins of the two electrons are opposite (Pauli paramagnetic 
limit). It is not clear what impact it has on Dirac semi-metals. If the impact is large it could be incorporated as an 
additional paramagnetic term in the GL energy. In an isotropic Dirac superconductor one has only one possible term 
in the GL energy term linear in paramagnetic coupling and consistent with symmetries: 

Fpar = NfipJ^iiA*xA)-B, (71) 

where is the effective ” spin” of the Gooper pair sometimes called ” Zeeman coupling” [l^ [ 2 ^ . The single particle 
Hamiltonian in magnetic field has the Pauli term • B, where the Bohr magneton, = e?i/2TOc, determines the 
strength of the coupling of the spin to magnetic field, with m being the free electron mass. The direct calculation, 
see Appendix A, shows that = 0. 

Acknowledgements. We are indebted to D. Li and M. Lewkowicz for valuable discussions. Work of B.R. was 
supported by NSG of R.O.C. Grants No. 98-2112-M-009-014-MY3 and MOE ATU program. 
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IX. APPENDIX A. LOCAL TERMS IN GL FREE ENERGY 
A. Gritical temperature calculation 

Starting from equation Ea. (l23ll the angle integrations result in (for /r >> Td,Tc) 

1 ^ _ 

3 np + + 2^2 {v%p^ + /r2) 

_ tanh {e/2T) ^ 2 Td7 ^ 

12n'^hPvp Jg=-TD ^ Gir'^hPvp ttTc 

where e = vpp — P- See the last (BCS) integral in[^. 


(72) 


B. Cubic terms coefficients calculation 

To fix the two coefficients, Pi and P 2 Ea. ipS)) we use only two components. The particular case j = k = I = 1 
(the coefhcient of gives after angle integration 


N{Pi+P 2 ) 


2T 


_ poo 

n “'P=0 


(A2) 


where 


S {p, n) 


Vpp^ {v%p‘^ + lOvpp^ {ujI - bp?) - 15 {p^ + 

+ 2n2,p2 (^2 _ ^2) (^2 ^2)2^ 


Performing finite integration (the upper bound on momentum, p + To, can be replaced by infinity), one obtains 


^ {Pi + P 2 ) 


8/^2 

Ibir'^T'^Vpt? 


(73) 


where the sum is 


S 3 


E 

n=0 


1 

{2n + lf 


7C(3) 

4 


Similarly taking j = l = 2,k = l (the coefhcient of A|2 Ai) gives after the angle integration 


NP 2 


2T 

157r2fi3 



Vpp^ (ryj-p"^ + 10u|,p2 3^2^ ^^2 

(vpP* + 2'u|.p2 (^2 _ ^2) (^2 ^2 )2^ 

4/i2 

Ib'n^v^pkPT'^ ^ 


(74) 


(75) 


resulting in Ea.dMl). 
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C. Effect of the Pauli interaction term 

The single particle Hamiltonian in magnetic field is 

K = —ivFhV-a. — /r + • B, (76) 

In order to fix the coefficient of the paramagnetic term linear in both the order parameter and Pauli coupling it 
is enough to expand the linearized Gorkov equations Eq. pij) to the first order in the spin density. Normal Greens 
functions have the following corrections: 


Dz ~ — HgD (E • B) D\ 

bz ^ 5 +(s‘• H) 5 . 

The Pauli term in Gor’kov equation (after multiplying by T/ and taking the trace as usual), Eq. 
obtained from expansion of Ea. ddTll . 


^Tr|r,DA*D| 


(77) 

, therefore is 

(78) 


Ljp 


Y2Z 

^ijk 


— i^,2;Eijk^j BkN — fj,QB^jf.AjBk, 

= ^ Tr {r,5 (eIDT* - T*D^k) d} , 


cjp 


The bubble sum is directly evaluated and vanishes = 0. 


X. APPENDIX B. CALCULATION OF GRADIENT TERMS IN THE GL 

Rotational invariance allows to represent the sum in Eq (1381) terms of coefficients ut and uf: 

— N {ut {P^Smj — PmPj) + ULPmPj) 

= PkPiY.Tr {t^^D (w, q) Tf U" (a;, q)} , 

a;q 

where 


D” = 


In particular 


and 


{qjai + qittj + Sij {iuj + fi + a ■ q) + 2qiqjD} . 


Nut. = - 


YTr{T,D{u;,p)T:D';Au^,P)} 


tup 






S3 


7C(3) /i2 7a^)ivFhfN 


607 r ^ T'^vf^i 


1207r2r2 


Nut = -Y P) P)} = 


(79) 


(80) 


(81) 


(82) 
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XI. APPENDIX C. COLLECTIVE MODES 


Here details of the calculation of free energy in harmonic approximation are given. The expansion of order parameter 
Eq.dm to quadratic order is 

^0 = (1 + ^) (^1; ^2; 1) = (0, 0, 1) + ((5l, ^2; s) + ^2^^ 0) . (83) 

The gradient terms that do not involve the vector potential A are: 


F. 


( 1 ) 


The terms involving A read 


Potential terms result in 


while the magnetic energy is: 


grad 

^0 


= Ut {djRadjRa + djiadjia + (djef'j 
I -Cdo 

+u ■ 


(84) 


+ 92 -^ 252-^2 + 92 / 292-72 + 


+2 


(92i?l) (9ii?2) + 

+ (92/i)(9i/2) + (9„e) (dsRc) 


^ = 2^^ (aMA„ + uAi, 


he 


F, 


pot 




= Ag|^l^e^-2/3,(/?+/|)| 


(85) 


( 86 ) 


A^ 


87re*2VA2^* A, 


(87) 


The fluctuation matrix M in Eg . (148)) thus was constructed. 


XII. APPENDIX D. EXACT SOLUTION EOR UPPER CRITICAL MAGNETIC FIELD 


In this Appendix the matrix R defined in Eg. (1561) determining the perpendicular upper critical field is diagonalized 


variationally. 


A. Creation and annihilation operators 

Using Landau creation and annihilation operators in units of magnetic length for the state with kx = 0 

(independent of x), so that covariant derivatives are 


Dx = dx + iy = iy = -^{a + a+) ; 

Dy = ■ 

In terms of these operators the matrix operator R takes a form: 

R = Ut + — + V; 

Vii = 2 uTa~^a + ^ (a^ + + 2 a~^a) ; 

Vi2 = V21 = — — a^) ; 

V22 = 2 uTa^a — ^ (a^ + — 2 a~^a) . 


( 88 ) 


( 89 ) 
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The exact lowest eigenvalue is a combination of two lowest Landau levels. Indeed applying the operator V on a 
general vector on the subspace gives 


/a|0)+/3|2) 

-(-te + ^+2/5)+4-r/3) |2) 

+ |(_/3_,^)|4) 

+ ^(*/3 + < 5 )| 0 ) 

V +fH/3 + <5)|4) y 


(90) 


For 6 = ip, higher Landau levels decuple and one gets eigenvalue equations 


—V —y/2u 0 

-4-T-2n-. 

0 uz-\/2 —V 


= 0 , 


(91) 


resulting in three eigenvalues of % 


Ml) 


ut + m/ 2, /i*-^) 

2>ut + —u ± 'sj 4My + Autu + ivP. 


(92) 
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Figures Captions 

Fig. 1. A complex vector field can be written as A = A (ncosx+tmsinx), where n and m are arbitrary unit 
vectors and 0 < x < 7r/2. 

Fig.2. Fluctuations of the order parameter in the unitary gauge can be parametrized by the five real fields 
A = Aq (1 + s) (i?i + ill, i ?2 + il 2 , 1) ■ 

Fig.3. 

Superconductivity arises from the normal state when the order parameter is formed in direction perpendicular to 
the magnetic field. 

Fig. 4 

The polarization vector of the incident beam rotates while passing a film of thickness d by angle (j) tan (j) = 

\B, {d)/By (d)|=exp[(A^i-Az')d]. 



